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ABSTRACT 

We examine the evolution of the Tully-Fisher relation (TFR) using a sample of 89 held 
spirals for which we have measured confident rotation velocities (V m t). This sample 
covers the redshift range 0.1 < z < 1, with a median of (z) = 0.33. The best-fitting 
TFR has a slope consistent with that measured locally, and we find no significant 
iy-j , evidence for a change with redshift, although our sample is not large enough to well 

constrain this. By plotting the residuals from the local TFR versus redshift, we find 
evidence that these luminous (Mb ^ MJj) spiral galaxies are increasingly offset from 
the local TFR with redshift, reaching a brightening of — 1.0±0.5 mag, for a given V TO t, 
£> ' by z ~ 1. This is supported by fitting the TFR to our data in several redshift bins, 

CN \ suggesting a corresponding brightening of the TFR intercept. Since selection effects 

would generally increase the fraction of intrinsically-bright galaxies at higher redshifts, 
we argue that the observed evolution is probably an upper limit. 

Previous studies have used an observed correlation between the TFR residuals 
and V I0 t to argue that low mass galaxies have evolved significantly more than those 
| with higher mass. However, we demonstrate that such a correlation may exist purely 

-»»« . due to an intrinsic coupling between the V rot scatter and TFR residuals, acting in 

' combination with the TFR scatter and restrictions on the magnitude range of the data, 

and therefore it does not necessarily indicate a physical difference in the evolution of 
q ■ galaxies with different V xo t- 

' Finally, if we interpret the luminosity evolution derived from the TFR as due to 

c/3 , the evolution of the star formation rate (SFR) in these luminous spiral galaxies, we 

find that SFR(z) oc (1 + z) 1 ' . Notwithstanding the relatively large uncertainty, 
this evolution, which is probably overestimated due to selection effects, seems to be 
slower than the one derived for the overall field galaxy population. This suggests that 
' the rapid evolution in the SFR density of the universe observed since z ~ 1 is not 

driven by the evolution of the SFR in individual bright spiral galaxies. 
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1 INTRODUCTION 

The sl ope, intercept and scat ter of the Tully-Fisher relation 
(TFR; iTulfv fc Fisherl Il977l) are key parameters that any 
successful prescription for galaxy formation and evolution 
must reproduce. As a relation between the luminosity and 
rotational velocity of disc galaxies, it is a combination of 
a number of fundamental galaxy properties and relations. 

* Based on observations made with ESO Telescopes at Paranal 
Observatory under programme IDs 066.A-0376 and 069.A-0312. 
f Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Institute. STScI is operated by the association of Uni- 
versities for Research in Astronomy, Inc. under the NASA con- 
tract NAS 5-26555. 
| E-mail: ppxspb@nottingham.ac.uk 



These include the relation between observed rotation ve- 
locity and total galaxy mass, that between the total and 
luminous mass distributions, and the mass-to-light ratio of 
the galaxy's stellar population. The first of these is depen- 
dent upon the degree to which the stars in the galaxy are 
rotationally supported and, along with the second relation, 
is determined by the dynamical galaxy formation process. 
The third component, the mass-to-light ratio in a particu- 
lar photometric band, is determined by the star formation 
history of the galaxy. 

The TFR therefore provides a wealth of information 
about the astrophysical processes involved in forming and 
maintaining spiral galaxies. In the past decade significant 
progress has been made toward understanding the origins 
of the TFR. This involves determining how the parame- 
ters involved conspire to maintain the relation, through self- 
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regulation of star- formation in the disk dSilklll99^) . while 
explaining the variations which lead to its intrinsic scatter. 
Modern simulations are able to reproduce the s lope and scat- 
ter of the /-band TFR Ce.g.. lKoda et aljl200oh . respectively 
identifying these with the natural range of mass and spin pa- 
rameters for dark matter haloes. However, reproducing the 
TFR intercept while matching other properties of the galaxy 
populatio n is currently b eyond the abilities of semi-analytic 
models (|Cole et alJfc'oOfl . 

As well as a goal to understand in its own right, the 
TFR is particularly useful as a benchmark with which to 
compare samples of galaxies, in order to examine the differ- 
ences between them. Using this method, we can gain insight 
into the evolution of disc galaxies by considering the varia- 
tion of the TFR with cosmic time. In this paper we present 
the B-band TFR of field spirals at redshifts 0.1 < z < 1, and 
examine its evolution over the past 8 Gyr. We then assess 
the implications of these results for the evolution of the star 
formation rate (SFR) in spiral galaxies. 

In Section 2 we briefly describe our observations, data 
reduction, and techniques for measuring the galaxy param- 
eters. We present the TFR of our field galaxy sample in 
Section 3, along with an examination of its evolution with 
redshift, and a comparison with other studies. This is fol- 
lowed by a discussion of the implications for SFR evolu- 
tion in Section 4, and finally we give our conclusions in Sec- 
tion 5. Throughout we assume the concordance cosmology, 
with tt A = 0.7, n m = 0.3 and Ho = TO km s- 1 Mpc- 1 . All 
magnitudes are in the Vega zero-point system. 



2 DATA 

The target selection, photometry and spectr oscopy are de- 
scribed thoroughly in iBamford et, al.l 1)20051 hereafter Pa- 
per I), and therefore only briefly summarized here. The pri- 
mary motivation for our spectroscopic observations was the 
comparison of spiral galaxies in five clusters at0.2<z<0.8 
with field spirals at comparable redshifts. This has been ex- 
amined in Paper I. In the process we have observed a large 
number of field galaxies, by themselves forming a useful sam- 
ple for examining the evolution of the more general spiral 
galaxy population. This is the subject of the current paper. 

2.1 Target selection 

For each field, galaxies were identified and their magni- 
tudes measured using SExtractor on our i?-band imag- 
ing. The galaxies to be observed spectroscopically were 
then selected by assigning priorities based upon the like- 
lihood of being able to measure a rotation curve. Higher 
priorities were assigned for each of the following: disky 
morphology, favourable inclination, known emission line 
spectrum, and available HST data. A priori knowledge 
of the galaxies' approximate spectral type varied depend- 
ing upon the availability of such data in the literature. 
The following sources were utilized for each fi eld, labelled 
by th e target cluster in that field: MS0440: iGioia et all 
Il998l: AC114: ICouch fc Sharpies! Il987t ICouch et alJll998t 
A370: IDressler et all Il999t ISmail et al] ll997l :CL0054: 
iDressler et all Il999l: ISmail et all Il997l: P -A Due private 
comm.; MS 1054: Ivan Dokkum et al.ll20ofl This element of 
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Figure 1. The redshift distribution of our target galaxies. The 
hatched histogram shows the distribution of all field galaxies ob- 
served with identifiable emission lines, and the filled area only 
those in our final TFR field sample. 

the selection is therefore highly heterogeneous, and certainly 
not ideal. However, while it makes an understanding of our 
selection function very difficult, this does reduce the possi- 
bility of any systematic effects. 

For each spectroscopic mask, slits were added in or- 
der of priority, and within each priority level in order from 
brightest to faintest 7?-band magnitude. The only reason 
for a particular galaxy not being included is a geometric 
constraint caused by a galaxy of higher priority level, or a 
brighter galaxy in the same priority level. Often the vast 
majority of the mask was filled with slits on galaxies in the 
highest priority levels, with occasional recourse to lower pri- 
ority objects in order to fill otherwise unoccupied gaps in the 
mask. The redshift distribution of our field sample is shown 
in Fig. □ 

2.2 Photometry 

The photometry for this study was principally measured 
on HST archive im ages, supplemented by our own FORS2 1 
JSeifert et alfeOOOD B-band imaging and additional ground- 
based data, kindly provided by Dr. Ian Smail. These lat- 
ter images provide additional colour information, improving 
the fc-correction by helping constrain the galaxy SED. The 
magnitudes used in this paper are measured w ithin SExtrac- 
tor AUTO (Kron-style) apertures l|Bertin fc Arnout Jl99fih . 
and are correc ted for Galactic extinction using the m aps and 
conversions of ISchleeel. Finkbeiner fc Davisl il998l) . 

The observed apparent magnitudes were con- 
verted to apparent rest-frame B-band magnitudes by 
a colour- and fc-correction, determined using the observed 
colour information whe r e ava ilable and the SEDs of 
lAraeon-Salamanca et all jl993T) . The absolute rest-frame 
B-band magnitudes were then obtained assuming the con- 
cordance cosmology. Finally the magnitudes were corrected 
for internal extinction (including face-o n extinction of 
0.27 m ag), following the prescription of iTullv fc Fouauei 
(1985), to give the corrected absolute rest-frame B-band 
magnitudes, Mb, used in the following analysis. Note that 
both the cosmology and internal extinction correction pre- 
scription were chosen to allow straightforward comparison 
with other recent studies. 

A plot of Mb versus redshift is shown in Fig. [5] No- 
tice how our sample is limited to brighter brighter mag- 
nitudes with increasing redshift. We can assess how much 
of the typical galaxy population we sample by comparing 
with the M* luminosity function parameter. iNorberg et al] 

1 http://www.eso.org/instruments/fors 

© 2005 RAS, MNRAS 000. HUTU 



The TFR of distant field galaxies 3 



1 1 1 1 1 1 1 


1 1 

) 




• + : 






/o 

1 





I _i I i I i I i I i L 

0.2 0.4 0.6 0.8 1 



z 

Figure 2. Absolute rest-frame B-band magnitude versus redshift 
for all our field galaxies with identifiable emission lines. Galaxies 
in our final field TFR sample are shown by filled points, while 
open points indicate those for which no emission lines pass our 
quality control criteria (see Section 12.31 . The dashed line indi- 
cates Mg as calculated in the text. The dotted line shows the 
absolute magnitude corresponding to an apparent rest-frame B- 
band magnitude of 22.5 mag, as a function of redshift in our 
adopted cosmology. 



(2002) use data from the 2dF Galaxy Redshift Survey to 
measure the 6j-band luminosity function, finding M 6 * 7 — 
51og(/i) = —19.66 ± 0.07 mag. We can tran sform this into 
the_B-band using the conversion suggested bv lNorberg et all 
( 2002), bj = B - 0.28(B - V), and the colour of a typical 
galaxy in our sample (estimated from the best-fitting SEDs), 
(B — V) = 0.52. With H = 70 km s" 1 Mpc" 1 , we therefore 
have M B ~ —20.3 mag. This is shown by the dashed line in 
Fig-IU At z ~ 0.2 we sample a range ~2 mag either side of 
M B . By 2 ~ 0.5 this has reduced to (M B - 2) < Mb < M% 
mag, and at z ~ 0.8 we are limited to Mb < (M B — 1) 
mag. We therefore sample most of the giant spiral popula- 
tion below z ~ 0.5, but beyond this we are limited to only 
the brightest galaxies in this class. Note that this effect is 
less serious than in conventional magnitude limited stud- 
ies as our galaxies at high redshifts have, on average, been 
observed with longer spectroscopic integrations. This is be- 
cause the majority of our high-redshift galaxies are from 
masks targeting our more distant clusters, and thus with 
longer exposure times (see Section [2.311 . 

Inclinations (i) and photometric scalelengths (ra lP hot) 
for the disk components of the observed galaxies were mea- 
sured in bands close to R, preferentially in the HST im- 
ages (b ands F606W, F675W or F702W), and primarily using 
GIm2d dSimard et all2002T) . Reliable HST inclination (scale- 
length) measurements were available for 47 (39) per cent of 
the galaxies; for the remainder inclination (scalelength) was 
measured on the ground-based i?-band imaging, again usu- 
ally by GIm2d. These measurements are therefore separated 
from the effect of the bulge component and corrected for the 
effect of seeing. For a small number of galaxies the GIm2d fit 
was unreliable. In these cases the SExtractor axial ratio 



was used, assuming an infinitely thin disk (as does GIm2d). 
These inclinations were corrected by a factor determined 
from an empirical comparison of SExtractor and GIm2d 
based inclinations. One factor was used for all ground-based 
measurements, and another for the HST measurements. 

In order to further evaluate our selection function, we 
measure effective radii, r e ff , of circular apertures containing 
half the galaxy light. These were obtained from our FORS2 
i?-band imaging, the only band which is available for nearly 
all of our galaxies, 2 using SExtractor's FLUX_RADIUS 
output. Plots of rd.phot and r e s versus redshift are shown in 
Fig. 13 Both of these plots include a dotted curve indicating 
the physical size, in kpc, corresponding to an observed an- 
gle of 1 arcsec, as a function of redshift. The seeing in the 
i?-band imaging is roughly 1 arcsec FWHM, but for compar- 
ison with r e f; one must convert this to a half-light radius. A 
Gaussian profile with 1 arcsec FWHM has a half-light radius 
of ~ 0.6 arcsec. The measured r e s are thus all larger than 
the seeing, but from the correlation of the lower bound of 
the points r e g distribution with redshift, the seeing is often 
dominating the measurement. 

In order to compare the seeing scale with rd lP hot one can 
note that if the exponential scalelength of a point source, 
with a Gaussian seeing profile, were measured, the result 
would be close to the radius at which a Gaussian profile 
reaches 1/e of its peak value. This occurs at a radius of 
0.6 times the FWHM. However, this comparison is compli- 
cated by the inclusion of a bulge component in the fit sur- 
face brightness model. To some extent, though, the fact that 
the rd, p hot distribution extends below the seeing scale, and 
shows little variation in its lower limit with redshift, indi- 
cates that rd.phot is rather more independent of the seeing. 

However, due to the similarity of the galaxy sizes to 
the seeing scale, and in particular the impact of limited 
resolution on the bulge-disc decomposition, our disc scale- 
lengths are potentially unreliable. This problem becomes 
significantly worse for both higher redshift and intrinsically 
smaller galaxies. We therefore choose not to base any infer- 
ences on our measurements of galaxy size, and only use them 
to give an indication of any biases in our sample selection. 

Another tool to examine our sample selection is pro- 
vided by the surface brightness of our galaxies. We calcu- 
late the apparent ii-band average surface brightness within 
r eS , which we denote p B , c ff, a pp, using /Uft, e s,a PP = -Reff + 
2.5 log 27rr 2 ff , where _R e ff — R + 2.5 log 2 is the magnitude 
within r c ff by the definition of the half-light radius, with to- 
tal magnitude, R. This, converted into an absolute observed- 
frame R-ba,nd surface brightness using the distance modulus 
for our adopted cosmology (denoted /iij,eff,abs), is plotted in 

Fig.m 

A comparison of Figs. |21 and |1] reveals the two plots to 
show very similar behaviour. This, along with the apparent 
lack of any redshift-dependent selection on rd lP hot demon- 
strated by Fig. EI a), and the simple explanation that the 
correlation in Fig. E^b) is due to seeing affecting the r e fi 

2 Note that, due to full two-band HST coverage, FORS2 
-R-band imaging was not used by the study of MS1054 by 
iMilvang-.Tensen et alJ J2003t) . We supplement our data set with 
these earlier observations, as described later, but do not have R- 
band effective radii or surface brightness measurements for these 
additional galaxies. 
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Figure 3. Plots against redshift of (a) photometric disc scale- 
length, r d spcc and (b) effective (half-light) radius, r ef j, both in 
kpc, for all our field galaxies with identifiable emission lines (ex- 
cluding the supplementary MS1054 field in panel (b)). Galaxies 
in our final field TFR sample are shown by filled points, while 
open points indicate those for which no emission lines pass our 
quality control criteria (see Section 12.31 . The dotted line indi- 
cates the physical size subtended by an angle of 1 arcsec, as a 
function of redshift in our adopted cosmology. The dashed line in 
panel (b) shows the same for an angle of 0.6 arcsec, approximately 
equal to both the half-light radius and exponential scalelength of 
a Gaussian with 1 arcsec FWHM. 



measurement, rather than any selection effect, implies that 
the surface brightness distribution is largely a result of a 
magnitude-limited selection, rather than any selection on 
galaxy size. 

The distributions in all these plots of the galaxies which 
make our final TFR sample, relative to those in the whole 
sample, are considered in the following section, after a de- 
scription of the final TFR sample selection criteria. 




Figure 4. Plot against redshift of the absolute observed-frame R- 
band surface brightness within r e g , for all our field galaxies with 
identifiable emission lines (excluding the supplementary MS1054 
field). Galaxies in our final field TFR sample are shown by filled 
points, while open points indicate those for which no emission 
lines pass our quality control criteria (see Section 12.31 . The dot- 
ted line indicates the absolute surface brightness corresponding 
to an observed apparent surface brightness of 25 mag arcsec -1 , 
according to our adopted cosmology. 



2.3 Spectroscopy 

The spectroscopy fo r this study was observed using FORS2 
JSeifert et al.ll200d) . in MXU mode, on the VLT. In this 
mode, multiple slits are cut into a mask which is placed in 
the focal plane. The slits were individually tilted to align 
with the major axis of each galaxy, as determined from the 
J?-band pre-imaging. Slits of 1 arcsec width in the disper- 
sion direction were used, with lengths that varied in order 
to sample sky at the edge of each galaxy, while efficiently 
packing the masks with targets. Exposure times varied from 
30 to 210 minutes per mask, depending on the redshift of the 
cluster that formed the primary target for each observation. 
The seeing was typically ~ 1 arcsec and always below 1.2 
arcsec. These data and their reduction are described more 
completely in Paper I, and are therefore only briefly sum- 
marised here. 

The data were reduced in the usual manner 
( Milv ang- Jensen! 120031 Bamford thesis in preparation) 
to produce straightened, flat-fielded, wavelength-calibrated 
and sky-subtracted 2d-spectra for each galaxy. From these, 
individual sky- and continuum-subtracted emission line 
'postage stamps' were extracted. The main emission lines 
observed were [OIIJA3727, H/3 and [OIII]AA4959, 5007, with 
Ha but no [OII]A3727 for nearby galaxies (z < 0.3, depend- 
ing upon the position of the slit in the mask). In order 
to measure the rotation velocity (Kot) and emission scale 
lengths (rd, sp ec) we fit each emission line independently us- 
ing a synthetic rotation curve m ethod based on elfit2d by 
ISimard fc PritchetJ Jl998lll999|) . and dubbed elfit2py. In 
this technique a model emission line is created for particular 
sets of parameters, assuming a form for the intrinsic rotation 
curve, an exponential surface-brightness profile, and given 
the galaxy inclination, seeing and instrumental profile. The 
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intrinsic rotation cur ve assumed here i s the ' universal rota- 
tion curve ' (URC) of iPersic fc Saluccil Jl99ll) . with a slope 
weakly parametrized by the absolute B-band magnitude, 
Mb- From a comparison of the results from fits using both 
flat and 'universal' intrinsic rotation curves, adopting a flat 
rotation curve leads to values of Vrot ~ 10 km s _1 lower than 
found with the URC. The choice of rotation curve does not 
appear to affect the conclusions of this study. A Metropo- 
lis algorithm is used to search the parameter space to find 
those which best fit the data, and to determine confidence 
intervals on these parameters. Images of model lines with 
the best-fitting parameters are also produced for compari- 
son with the data. 

In order to produce a single value of V TO t and rd, spe c for 
each galaxy, the values for the individual lines are combined 
by a weighted mean, as described in Paper I. 

A significant fraction of the emission lines identified 
display dominant nuclear emission, or asymmetries in in- 
tensity, spatial extent or kinematics. In severe cases these 
departures from the assumed surface brightness profile and 
intrinsic rotation curve mean that the best-fitting model is 
not a true good fit to the data. A similar situation can occur 
for very low signal-to-noise (S/N) lines, where an artifact 
of the noise overly influences the fit. More concerning is the 
case of very compact lines, where the number of pixels is 
on the order of the number of degrees of freedom in the 
model, and hence an apparently good fit is obtained despite 
a potentially substantial departure from the assumed surface 
brightness profile. 

In order to eliminate such 'bad' fits a number of qual- 
ity tests are imposed, based on a measure of the median 
S/N and a robust x 2 goodness-of-fit estimate. In addition, 
the emission lines were 'traced' by fitting a Gaussian to each 
spatial column of pixels, and the region determined for which 
the trace is reliable. The distance from the continuum cen- 
tre to where the line could no longer be reliably detected 
above the noise we term the extent. Additionally, quantities 
describing the asymmetry, in terms of extent and kinemat- 
ics, and the flatness of the line at maximum extent were 
formulated. Unfortunately, a satisfactory set of quantitative 
criteria, to exclude 'bad' fits while efficiently retaining 'good' 
fits, could not be found. Therefore, we supplemented cuts on 
S/N and \ 2 with a visual inspection of every emission line, 
supported by the various quality tests just described. In the 
process galaxies with apparently dominant nuclear emission 
were excluded, as were lines producing zero Vrot estimates 
and those clearly inconsistent with a number of other lines 
for the same galaxy. 

As in Paper I, we have supple mented our data with that 
from lMilvang- Jensen et"aH fcOOSlh This is also FORS2 MXU 
data, reduced in a very similar way to that described above. 
In order to consistently combine the MS 1054 data with this 
sample, the emission line postage stamps for the MS1054 
galaxies have been re-fit using elfit2py and the same line 
quality criteria applied. 

After all these quality checks, our full field galaxy TFR 
sample contains 89 galaxies, with a mean of ~ 2.3 emission 
lines contributing to the measurements for each galaxy. This 
sample covers the redshift range 0.1 < z < 1, with a median 
redshift of (z) = 0.33, as show in Fig. 

Figures [5] [3] and 2] show the distributions of Mb, 
fd.phot, ''eft and /ifl, c ff,abs versus redshift for all field emis- 




log(V rot [km/s]) 

Figure 5. The Tully-Fisher rel ation for our full field TFR sam- 
ple. The fiducial local relation of lPierce &: Tullvl <1992l) is marked 
by the dot-dashed line, with its 3<r intrinsic scatter delimited by 
dotted lines. A weighted least-squares fit to all the points is shown 
by the solid line, see the text for details. The error bars in the top 
left corner indicate the 10th-, 50th- and 90th-percentiles of the 
(broadly unimodal) distribution of uncertainties on the displayed 
points. 

sion line galaxies observed. Those in our final TFR sample 
are marked by filled points, while those which are not, i.e., 
for which none of the emission line fits passed our quality 
criteria, are shown by open points. There is no noticeable 
difference between the distributions of selected and rejected 
galaxies in these plots. This indicates that our quality se- 
lection criteria are not biasing our sample from the point 
of view of the galaxies' broadband photometric properties, 
beyond those baises inherent to our initial spectroscopy ob- 
servational selection procedure. 



3 THE TULLY-FISHER RELATION 
3.1 Basic fit 

Our TFR for distant field galaxies is given in Fig. The 
solid line is an 'inverse' fit to the data, 

log Vrot = a + bM B (1) 

found by minimising the weighted squared deviations in 
log Kot • The weights applied to each point are 

t«i = l/o? , (2) 
where 

&i = ofog Vtot: i + b 2 al lBti + of nt . (3) 

Here oi og y rot ,i and om b ,i are the errors on each point derived 
from the reduction process, and o"i n t is the intrinsic scatter of 
the TFR, which is allowed to vary such that the reduced chi- 
squared statistic, x?, is unity. This was achieved by iteration, 
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each time recalculating the weights using the new value of 
b, and <7i n t determined via the recurrence relation 



(4) 



where, as usual for two degrees of freedom, the reduced-chi- 
squared statistic is 



Xr 



X' 



v-^ / l°g Vroti — a — bMs.i 



n-2' ^" = Z.l- — J • ( 5 ) 



The (weighted) total scatter, <T tota i, is calculated by 



2 

""total 



^2 Wi o> 



Og Vrr 



bM B ,i 



(6) 



Notice that when x r is a t the desired value of unity, the 
recurrence relation (Eqn.[IJ reduces to of ntJ+1 = crf nt 
as we would want. When Xr > 1 this implies of nt is too small, 
and thus it is increased for the next iteration. Alternatively, 
when x? < 1 this implies af nt is too large, and it is therefore 
decreased for the next iteration. The value of a may be tuned 
to minimise the number of iterations before convergence is 
achieved; a = | was found to work well. 

The fit converges in ~ 4 iterations to give the TFR 
(converted back into the 'forward' form): 



M B = (-8.1 ± 0.8) • log Vrot + (-3.2 ± 2.0) 



(7) 



with ototai = 1.0 mag and c^nt = 0.9 mag. 

The 'inverse' fit is often used to avoid a potential bias 
because of the apparent magnitude limited nature of a TFR 
sample. At faint apparent magnitudes, points are preferen- 
tially selected above the real TFR. For samples over a nar- 
row range of distance moduli, e.g., in local cluster studies, 
the nature of this bias is to cause an apparent flattening of 
the TFR slope. While th e depe ndence of V ro t on Mb still 
causes an effect, IWillickl il994h has shown that using an 
'inverse' fit reduces the bias by a factor of five compared 
with the conventional 'forward' fit. As this study covers 
a wide range of distance moduli, apparently faint galaxies 
are found over the whole range of absolute magnitude, de- 
pending on their redshift. Any bias on the best-fitting TFR 
should therefore be weakened, and take the form of an inter- 
cept shift rather than a flattening of the slope. This is par- 
ticularly true given the more random priority-based, rather 
than magnitude-based, selection. However, to minimise any 
such effect, and allow us to confidently use the same fitting 
method to compare differently selected samples, we choose 
to work with the 'inverse' fit. 

We weight our fit to make full use of the data, and 
avoid the influence of unreliable points, while the inclusion 
of an intrinsic scatter term prevents points with small errors 
from dominating the fit, and allows us to estimate this useful 
parameter. 

The thin lines in Fig. |K| indicate a fid ucial local field 
TFR This is derived from the TFR of iPierce fc TuUvl 
il992L hereafter PT92), with a zero-point adjustment be- 
cause PT92, while otherwise using the internal extinction 
correction of lTullv fc Fouauel(|l985l i. do not include the 0.27 
mag of face-on extinction that is applied to our data. The 
fiducial PT92 TFR, adapted to our internal extinction cor- 
rection, is thus: 



Ml T92 (K c 



-7.48 losKc 



3.37 



(8) 



PT92 obtained this relation by an 'inverse' least-squares 
fit, minimising the residuals in rotation velocity. This is the 
same method we have used to fit our TFR, except that our 
fit is weighted and includes the intrinsic scatter in these 
weights. The weighting should not cause any bias, and hence 
our fitting method is comparable to that used to produce the 
fiducial TFR of PT92. 

There is clearly a significant offset between the TFRs 
of PT92 and this study. This may not be entirely an evolu- 
tionary effect, and some part of it is likely due to the differ- 
ent manner in which we measure the rotation velocities and 
magnitudes compared with PT92. Note that in generating 
the fiducial relation we have assumed that V ro t = 0.5Wji, 
where W l R is the fully-corrected HI velocity width measured 
by PT92. 

A further issue is whether or not the absolute calibra- 
tion of the PT92 TFR is correct. This is based upon Cepheid 
and RR Lyrae distances to six local calibrator galaxies, and 
hence dependent upon the calibration of the whole extra- 
galactic distance scale at the time. There is also the issue 
of how representative this small number of local calibrator 
galaxies are compared with the whole PT92 TFR sample. 
iPiercd (Il994l . hereafter P94) uses the PT92 TFR to calibrate 
supernovae distances and derive Ho = 86±7 km s _1 Mpc -1 . 
If, as is suggested, the primary cause of uncertainty in this 
value is the absolute TFR calibration, then we could poten- 
tially use the modern, significantly more accurate value of 
Ho to correct this. We could therefore derive a more accu- 
rate TFR absolute calibration by using our cosmologically 
measured Ho, a reversal of the traditional method. 

The current best estimate of Ho is provided by a combi- 
nation of data from WMAP and a number of o ther surveys: 
H = 71+i km s- 1 Mpc- 1 dSpergel et al.ll2003i) . The differ- 
ence in distance modulus between this Ho and that derived 
by P94 is 0.5 ±0.1 mag, in the sense of moving the TFR to 
brighter magnitudes. The errors given only include those on 
the WMAP Ho measurement, as the P94 errors are dom- 
inated by the uncertainty on the distance scale, which we 
are replacing. It can therefore be argued that a correction 
of —0.5 mag should be applied to the intercept of the PT92 
TFR. 

Fortunately in this study the uncertainty on the inter- 
cept of the fiducial local TFR is of little concern. Our sample 
covers a wide range in redshift, selected and analysed in a 
homogeneous manner. We can therefore examine the TFR 
evolution using only this study's data, without recourse to 
external work. 

In future studies it may b e wise to consider usin g the 
more recent B-band TFR of iTullv fc Pierce! §000) as a 
comparison. This relation has an absolute calibration based 
on Cepheid distances to 24 galaxies and implies a value of 
H = 7 7 ± 8, more consi stent with the WMAP result. The 
work of IVerheiienl |2001) provides another more recent lo- 
cal comparison TFR, which has already been used by some 
groups. However, th ese studies use th e internal extinction 
correction scheme of lTullv et alJ il99Sf) . which has a strong 
dependency upon galaxy luminosity (or alternatively Vrot)- 
While this form of the internal extinction should be more 
accurate, it will require care as the dependency on luminos- 
ity (or Vmt) is calibrated locally, and may not be valid for 
distant galaxies. 
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The internal scatter we measure, <7i n t = 0.90 mag, is 
considerably larger than the ~ 0.4 mag generally found fo r 
local samples (e.g iPierce fc Tullvl 11993: iDale et ai"lll99Sl) . 
However, most local studies are focused upon using the 
TFR as a distance indicator. They therefore impose very 
strict selection criteria, e.g., requiring undisturbed late- type 
spirals, in an effort to produce as strong a correlation as 
possible. At high redshift we do not have the luxury of 
abundant data, and so cannot impose such strict criteria. 
For example, it is harder to identify galaxies with slightly 
disturbed rotation curves and morphologies, making rejec- 
tion of such objects impossible. In addition, corrections for 
known, locally-calibrated correlations, e.g., between TFR in- 
tercept and morphological type, are often applied to local 
samples. However, such corrections are not yet calibrated 
at high-redshift and so cannot be applied to our data. The 
variations caused by minor disturbances, differing sample se- 
lection, and varying corrections may be responsible for the 
entire increase in intrinsic scatter. However, physical effects, 
such as an increased stochasticity of star-formation, may 
also contribute. 



3.2 TFR evolution with redshift 

The slope of the TFR we measure for the whole field galaxy 
sample is very similar to that found locally by PT92, and 
consistent within our errors. Note that we have used a com- 
parable fitting method to PT92, so this statement is valid. 
The precise slope measured is very much dependent on the 
fitting method employed, and may explain why some studies 
find apparently different slopes. Our fit suggests there is lit- 
tle change in the TFR slope with redshift. We can attempt 
to evaluate this further by fitting redshift sub-samples of our 
data. 

As noted above, any difference in intercept between our 
sample and the local relation of PT92 may not be a real 
effect. However, by fitting redshift sub-samples we may in- 
vestigate any evolution of this offset purely within our own 
sample. 

In order to examine the intercept evolution we divide 
the sample in to five bins, each of 18 points (17 for the high- 
est redshift bin). The TFR is fit for each bin using the usual 
method, but with the slope constrained to the PT92 value 
(and hence very similar to the value for our whole sample) 
to isolate the changes in intercept. In panel (a) of Fig.|S]wc 
plot these results, as the Mb value of the best-fitting TFR 
at the median Kot of the full sample. Also plotted is a con- 
ventional weighted least-squares fit to the points. There is 
reasonably clear evidence for a brightening of the TFR in- 
tercept with redshift of ~1 mag by z = 1, corresponding to 
a factor of ~2.5 in luminosity at a fixed rotation velocity. 

At high redshift, however, we only sample the bright end 
of the galaxy luminosity function, and hence preferentially 
select objects which have been brightened. In particular we 
require a certain emission line flux to measure Kot, and 
hence the limiting emission line luminosity of our sample will 
increase with redshift. As emission line luminosity increases 
with SFR, we will therefore sample galaxies with higher SFR 
at higher redshifts, which in turn implies brighter B-band 
magnitudes for the high-redshift galaxies. The evolution of 
the TFR intercept that we measure is therefore probably an 
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Figure 6. Plots showing the evolution of the slope and intercept 
of the best-fitting TFR in several redshift bins, indicated by the 
horizontal error bars. Panel (a) gives the intercept (at the median 
Vrot of the full sample) for TFR fits to five redshift subsamples of 
the data. Each redshift bin contains ~ 18 points. The line shows 
a weighted least-squares fit to the points. For these fits the slope 
was held fixed to the PT92 value. Panel (b) plots the slope for 
TFR fits to three redshift subsamples of the data. Each redshift 
bin contains ~ 30 points. The horizontal lines show the local TFR 
slope of PT92 (dotted) and the slope from Eqn. [7] the fit to our 
whole field TFR sample (dashed). 



upper limit on the true brightening, e.g., at a given galaxy 
mass. 

To investigate any evolution in the TFR slope we divide 
the full sample into only three redshift bins, as the slope 
requires more points to constrain it. Each bin thus contains 
30 points (29 for the highest redshift bin). A TFR was fit to 
each sub-sample in the same way as for the full sample, and 
the results are shown in panel (b) of Fig. |S| Due to the small 
numbers of points in each bin, the errors on the slope are 
quite large. Given these errors, and the restricted Mb range 
of the data in the highest redshift bin, no strong constraints 
can be inferred. However, the slope in each bin is consistent 
with no evolution of this parameter with redshift. 

In order to further investigate changes in the intercept 
of the TFR via the offsets of individual galaxies, it is helpful 
to work with the residuals from the fiducial TFR: 



Mb - Ml™ (l/ rot ) 
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Figure 7. The residuals, AM B TF , from the fiducial PT92 
TFR plotted against rcdshift and look-back time. Galaxies with 
Mi? eff abs < —20.0 are denoted by squares, those with —20.0 < 
Mi? eft abs < —18.0 are shown by triangles, circles denote those 
with /i R c ff abs > —18.0, and galaxies with no ii R e g at>s measure- 
ment (i.e., those from the supplementary MS1054 data) are indi- 
cated by diamonds. The fiducial local relation of PT92 is again 
marked by the thin dot-dashed and dotted lines. A weighted least- 
squares fit to all the points is shown by the solid line, see the text 
for details. The error bars in the top left corner indicate the 10th-, 
50th- and 90th-percentiles of the (broadly unimodal) distribution 
of uncertainties on the displayed points. 

= Mb - (-7.48 log V ro t - 3.37) . (9) 

We can evaluate an evolution of the TFR with redshift by 
looking at a plot of these residuals versus redshift, as shown 
in Fig. A trend in AMb with redshift, such that more 
distant galaxies are brighter for a given rotation velocity, 
is apparent. Fitting the data in a similar manner to that 
described for the TFR in Section 13.11 (minimising residuals 
in AMg TF ) produces the relation 

AM B TF = (-1.0 ± 0.5) ■ z + (0.8 ± 0.2) mag (10) 

with Ototai = 0.9 mag and cr- m t = 0.8 mag. This is in good 
agreement with the more qualitative result from Fig.|S] Once 
again, however, this is most likely an upper limit on the true 
luminosity evolution experienced by giant spirals in the field, 
due to our preferential selection of the brightest objects at 
high redshift. 

In Fig. |7| there is clearly some asymmetric scatter to 
brighter offsets, which decreases with redshift. This may be 
due, at least in part, to the fact that there is a larger scat- 
ter in the TFR around V TO t ~ 100 km s _1 . Only brighter 
galaxies, and hence those with V to t > 150 km s _1 and thus 
lower TFR scatter, are observed at high redshifts. However, 
this does not seem able to account for the whole effect, and 
there is a hint that this reduction in scatter with redshift is 
real. 

The points in Fig. |7| are shown using different sym- 



bols depending upon absolute surface brightness, /^eff.abs- 
There appears to be no obvious difference in the AMb TF 
distribution of high and low surface brightness galaxies at 
a given redshift. This argues against the variation in scat- 
ter being simply a consequence of sampling different surface 
brightness ranges at different redshifts. It also implies that 
the varying surface brightness selection function has rela- 
tively small impact on the evolution we measure. 

As mentioned in the previous section, there is a system- 
atic offset between our TFR and that of PT92, even when 
the evolution with redshift is taken into account (as indi- 
cated by the zero-point of Eqn. 1101 being inconsistent with 
zero). The earlier discussion of the uncertainties in the in- 
tercept of the PT92 TFR also applies here. Note that the 
'cosmological' correction of —0.5 mag suggested in Section 
l3.1l would bring the zero-points into much closer agreement. 
In addition, there are a variety of other possible explana- 
tions for the offset. For example, the asymmetric scatter to 
brighter Mb for a given V ro t , seen particularly at lower red- 
shifts, may influence the fit away from the PT92 relation. 
This scatter may be due to galaxies that would not have 
been included in the PT92 sample, for example due to in- 
teractions, or that are less prevalent locally. A more rapid 
evolution in the past ~2 Gyr, which is not well constrained 
by our data, could also be responsible. 

3.3 TFR residuals versus V TO t 

It may be thought that plotting the TFR residuals, AMb TF , 
versus Vrot would provide a test for a change of slope. How- 
ever, this must be treated very carefully, as the two vari- 
ables are intrinsica lly correlated through the scatter in V TO t- 
iBohm et alJ (|2004l . hereafter B04) use this plot to argue that 
galaxies with lower Vrot, and hence lower mass, are offset fur- 
ther from the local TFR than more massive galaxies. How- 
ever, a correlation between AMb TF and Vrot does not neces- 
sarily imply such a result. The fact that the two variables are 
not independent (cf. Eqn. [5J , combined with scatter in the 
TFR and restrictions on the magnitude range of the data, 
causes an intrinsic correlation, even in the absence of any 
true difference in TFR slope. 

In order to show this in more detail, Fig. |S] gives a 
pedagogical illustration of the effect of scatter and mag- 
nitude cuts on the distribution of points in the TFR and 
AMb TF - log Vrot plots. Panels (a) and (b) show three cor- 
responding lines in the TFR and AM_g TF - log Vrot plot, re- 
spectively. The thick solid line is an example ideal TFR, 
corresponding to the PT92 relation with a constant magni- 
tude offset of —1.06 mag, the median AMs TF of our data, 
in order to allow more direct comparison of the examples 
with our data. This example TFR is limited to the range 
1.9 < log V ro t < 2.3, and hence -18.64 > M B > -21.63. 
The thick dashed and dotted lines indicate loci of constant 
Mb at the faint and bright limits of the example TFR, re- 
spectively. In panel (b), the plot of AMs TF versus log Vrot, 
the example TFR becomes a (relatively short) horizontal 
line, while the loci of constant Mb are lines with a slope 
negative that of the example TFR, as is obvious from con- 
sideration of Eqn. |5| 

To show the transformation of distributions of points 
in the TFR plot to the AMs TF - log V ro t plane we construct 
samples of galaxies with simulated properties. These are de- 
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Figure 8. A simple demonstration of the intrinsic correlation between the TFR residuals from the fiducial PT92 TFR, ie., AMj , and 
Vrot- The fiducial local relation of PT92 is again marked by the thin dot-dashed and dotted lines. Panels (a) and (b) show instructive 
diagrams of the TFR and AM B TF versus log Vrot plot, respectively. Each line in panel (b) corresponds to the same style of line in panel 
(a). Panels (c) and (d) show the distribution of points, as greyscale and contours, for simulated galaxies with true log(V _ ro t) uniformly 
distributed on the range 1.9-2.3, true Mg assigned to each galaxy corresponding to the PT92 TFR, with a constant offset to make it 
more easily comparable to our intermediate redshift data. Gaussian scatter with 0.133 dex standard deviation (corresponding to 1.0 mag 
in terms of Mg assuming the PT92 TFR slope) has been added to generate the observed log V ro t- No scatter is added to the Mg values. 
Panels (e) and (f) show the distribution of simulated galaxies with true properties created in the same manner as the previous panels. 
However, in generating the observed values in this case, no scatter has been added to log Vrot, while Gaussian scatter with 1.0 mag 
standard deviation is added to the Mg. In addition, to demonstrate the effect of a magnitude cut, points with observed Mg > —18.64 
have been removed from the sample. See the text for a more detailed discussion. 



signed to cover parameter ranges fairly similar to our data, 
but in order to keep our arguments straightforward, we ini- 
tially keep these simulations highly simplified. More realistic 
simulations are considered later. 

Each simulated galaxy is assigned a 'true' rotation ve- 
locity from a random uniform distribution in log Kot be- 
tween 1.9 and 2.3, and a 'true' magnitude corresponding 
to this 'true' rotation velocity using the example TFR de- 
scribed above. The 'true' galaxy properties are thus uni- 
formly randomly distributed on the solid line in panel (a) of 
Fig. |H| The distribution of 'true' properties is the same irre- 
spective of whether the initial property assigned is log Kot 
or Mg. Scatter, in terms of measurement errors or intrin- 
sic scatter in the TFR, may then be added to either the 
'true' rotation velocities, magnitudes, or both, to produce 
simulated 'observed' values of the parameters. 

The effects of measurement error and intrinsic scatter 
on the distribution of points are equivalent for these exam- 
ples. Differences between the two can only arise as artifacts 
of the simulation method. Measurement error must always 



be applied after both parameters have been assigned, while 
intrinsic scatter may be applied between the production of 
the first parameter, e.g., log Kot and the assignment of the 
corresponding second parameter, e.g., Mg. If the intrinsic 
scatter is only applied to the value of the first-assigned pa- 
rameter that is then used to generate the second parameter, 
then this is equivalent to a measurement scatter applied di- 
rectly to the second parameter. If the intrinsic scatter is only 
applied to the first parameter after the second parameter has 
been generated, then this is equivalent to measurement error 
on the first parameter. Alternatively, if the intrinsic scatter 
is applied to the 'true' value of the first-assigned parame- 
ter, which is then used to generate the second parameter 
and also used as a basis for the 'observed' value of the first 
parameter, then this is effectively an alteration of the ini- 
tial distribution from which the first parameter is derived, 
as well as the addition of a scatter in the second param- 
eter which is indistinguishable from a measurement error. 
As we wish to be able to control the distribution of the 
'true' parameters independently of the inclusion of scatter, 
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Figure 9. A more detailed demonstration of the intrinsic correlation between AMb TF and Vrot. The fiducial local relation of PT92 
is again marked by the thin dot-dashed and dotted lines. Panels (a) and (b) show the distribution of points in the TFR and AAfg TF 
versus log V TO t plot, respectively, for a realistic simulated sample of galaxies, as described in the text. This may be compared with panel 
(c), a plot of AMj TF versus log Kot for our real field TFR sample. Note the clear correlation, which, as shown by the simulations, may 
be explained purely by the intrinsic correlation between AMb TF and log Vrot , along with the effect of scatter in the TFR combined with 
a restriction on the magnitude range of the data. This correlation, therefore, does not necessarily imply that the data prefers a slope of 
the TFR different to that of the chosen fiducial TFR, and is thus not adequate evidence for such a real, physical effect. 



we choose to implement intrinsic scatter as offsets applied 
following the generation of both 'true' parameters, in which 
case it is identical to a scatter due to measurement error. 
For the following simple examples we therefore do not dis- 
tinguish between the two. All random scatters applied in 
these simulations are based on Gaussian distributions with 
the specified standard deviation. 

In panel (c) of Fig. |H1 we plot the distribution of 
10 s simulated points, as greyscale and contours, generated 
by the procedure described above, and with the inclusion 
of 0.133 dex of scatter in log Kot (corresponding to 1.0 
mag in terms of Mb assuming the PT92 TFR slope). As 
no scatter has been added to the magnitude coordinates, 
the distribution is restricted to a sharply defined range of 
Mb- Panel (d) shows the corresponding distribution in the 
Ail/g TF - log Vrot plane. The extension of the points paral- 
lel to the fiducial TFR produces a horizontal extension in 
this plot. The scatter in log Kot spreads the distribution in 
both the horizontal and vertical direction, along lines with 
slope equal to minus the TFR slope. The restrictions in Mb, 
due to the initial choice of the 'true' TFR points' distribu- 
tion, translate into sharply defined sloping cut-offs in the 
AMb TF - log Vrot plot, corresponding to the dashed and dot- 
ted lines shown in panels (a) and (b) of the same figure. 

Panel (e) of Fig. |S] shows, in the same manner as panel 
(c), the distribution of another 10 5 simulated points. These 
points were generated by the same procedure, but this time 
with 1.0 mag of scatter in Mb, and none in log Vrot. In addi- 
tion, in this panel we demonstrate the effect of a magnitude 
cut, by rejecting points with observed Mb > —18.64 (cor- 
responding to the faint end of the simulated 'true' Mb dis- 
tribution) from the sample. Without a magnitude cut, the 
edges of the distribution are defined by sharp vertical lines 
of constant log Vrot, due to the initial range of the simulated 
distribution, and by softer drop-offs parallel to the example 
TFR, due to the introduced scatter. 

In the AMs TF - log Kot plot, shown in panel (f), this 
translates into a rectangle, the constant log Vrot edges are 
obviously preserved, and the scatter is purely in the AAfs TF 



direction. However, the inclusion of a maganitude cut causes 
a slope at the corresponding edge of the AA/g TF - log Vrot 
distribution. For a more realistic, magnitude-limited sample, 
the range in log Kot would not be restricted as it is in this 
example, and so the entire low log Vrot side of the AMb TF 
versus log Kot distribution would be sloped. In addition, a 
realistic magnitude distribution would show a cut-off at the 
bright end due to the observed nature of the luminosity 
function. Magnitude errors in studies like this one do not 
typically exceed 0.2 mag, and so adopting 1.0 mag of Mb 
scatter would require intrinsic scatter in the TFR. However, 
assuming a uniform distribution of log Kot with this amount 
of intrinsic scatter is clearly at odds with the observed lu- 
minosity function. 

With an unrestricted range of log Kot and magnitude 
cuts for Mb > -18.64 and Mb < -21.63 the distribution 
of points in the TFR and AM b tf - log Kot planes would 
look nearly identical to those shown in panels (c) and (d). 
Therefore, scatter in either Mb and/or log Kot, in the form 
of intrinsic scatter or measurement errors, when combined 
with restrictions on the Mb range, whether due to the un- 
derlying luminosity function, selection effects or directly im- 
posed, produces a correlation in a plot of AMb TF versus 
log Kot • As these restrictions only affect the edges of the 
AMb TF - log Kot distribution, the strength of the correla- 
tion, and its apparent slope, are related to the ratio of the 
total scatter to the Mb range covered. The correlation is 
strongest, with slope equal to minus that of the TFR, when 
the Mb range is smaller than the scatter, and weakens, with 
flatter slope, as the Mb range grows relative to the scatter. 

Finally, note that the combination of a scatter with re- 
strictions on the magnitude range leads to the measured 
slope of the simulated/observed TFR appearing flatter than 
the real underlying TFR, as can be seen in panel (c) of Fig.|H| 
(Likewise, restrictions in the rotation velocity range lead to a 
steepening of the apparent TFR.) However, this well known 
problem is greatly reduced by choosing a method to fit the 
TFR which minimises the residuals in log Kot , rather than 
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M B (e.g., see lWillicklll99i) . as is done throughout this pa- 
per. 

While we have made the reasonable assumption of 
Gaussian scatter on log Kot and Mb in the above exam- 
ples, the appearance of a correlation does not depend on 
the form of the scatter, merely that there is some broad- 
ening of the TFR. Clearly the detailed distribution in the 
Ail/s TF - log Vtot plane will vary with the choice of the scat- 
ter distribution, but the existence of a significant correlation 
is robust. 

The examples in Fig. |H1 are relatively simple, in order 
to avoid confusing the effects of different simulation param- 
eters. However, to compare more fairly, though still quali- 
tatively, with our real data, we require a more detailed and 
realistic set of simulated data. To generate this data we use 
the following method. We assign each simulated galaxy a 
redshift, randomly drawn from a uniform distribution 3 on 
the range 0.0-1.0. Because of the varying volume elements 
enclosed by a given solid angle at different redshifts, we must 
preferentially include galaxies in our sample at a rate pro- 
portional to the comoving volume element at their redshift. 
To do this we calculate the ratio of the comoving volume 
element at the assigned galaxy redshift to that at redshift 
1.0. If a random number between and 1 is greater than this 
ratio, then we reject this galaxy from our sample and begin 
generating a new galaxy. This assumes that the comoving 
number density of galaxies is constant, which is true in the 
absence of mergers, and a reasonably good assumption at 
these redshifts. Changes in the number density at a given 
magnitude due to luminosity evolution are considered later. 

If the galaxy is not rejected it is then assigned a true 
absolute magnitude, randomly drawn from a uniform dis- 
tribution between —12.0 and —24.0. We model the effect of 
the luminosity function by calculating ^>(Mb), normalised 
to be less than unity on the magnitude range considered, 
and rejecting the galaxy from our sample if a random num- 
ber between and 1 is greater than 4>{Mb)- The l uminosity 
function assumed, 4>{Mb), is a Schecter function dSchechterl 
I1976T) with parameters from th e 2dF Galaxy Redshi ft Sur- 
vey (2dFGRS) determined bv iNorberg et alJ J2002I) . with 
M* converted from the fej-band to the _B-band as described 
in Section l2~2l i.e., M* = -20.3 and a = -1.2. 

An added complication for comparing with the real 
data is evolution of the luminosity function. Without in- 
cluding any evolution, the bright end of the luminosity func- 
tion of the simulated galaxies does not match that seen for 
our real sample. Although our data is not a complete sam- 
ple, and is subject to complex selection effects, the bright 
magnitude cut-off should be reasonably close to that of the 
true underlying luminosity function. In Section [3.21 we find 
evidence for luminosity evolution, at a given rot ation ve- 
locity, of u pto around — 1 mag p er unit redshift. lLovedavl 
(2004) and iBlanton et alJ i2003ft find evolution in the r- 
band luminosity function from SDSS data, although with 
only z < 0.3, amounting to a change in M* between —1 
and —2 mag per unit redshift. Evolution of the luminosity 
function with redshift is therefore certainly plausible, and 



3 actually discrete values separated by 0.001 to improve the com- 
putational speed, by allowing redshift-dependent values to be cal- 
culated prior to the galaxy-generating loop 



would have most noticeable effect around the sharp cutoff 
at bright magnitudes. To account for this we apply an evo- 
lution of M*(z) = M*(0) - 1.0 z when calculating tj>(M) for 
each simulated galaxy. This brings the luminosity function 
of the simulated galaxies into much better agreement with 
that of our real data. 

For the galaxies remaining, which now have the red- 
shift distribution expected given the comoving volume of 
space sampled at each redshift, and an absolute magnitude 
distribution corresponding to a realistic, evolving luminosity 
function, we now generate apparent magnitudes, and apply 
a magnitude limit. 

No intrinsic scatter is added to the magnitudes as we 
have already used a realistic observed luminosity function to 
generate the magnitudes. Adding intrinsic scatter to these 
values does not make sense, as it would change the assumed 
magnitude distribution to be less representative of that ob- 
served. Any intrinsic scatter in the TFR must therefore be 
added to the rotation velocities. 

To apply a realistic magnitude limit we first calculate 
'true' apparent magnitudes from the absolute magnitudes 
using the distance modulus at each galaxy's redshift. We 
then reject any galaxy with a 'true' apparent magnitude 
fainter than a specified limit. We apply a magnitude limit 
similar to that seen in our data, from examination of the plot 
of Mb versus z in Fig. [5] and comparison of the faint end 
of the luminosity functions for our simulated and real data. 
An apparent magnitude limit of 22.5 mag is thus adopted. 

Next, simulated 'observed' Mb are produced, by adding 
to the 'true' magnitudes values randomly drawn from a dis- 
tribution with properties similar to the measurement uncer- 
tainties on our real data. Rather than overcomplicate the 
simulation by allowing the measurement uncertainites to 
vary from galaxy to galaxy, we model the magnitude mea- 
surement scatter by a Gaussian distribution with a constant 
standard deviation of 0.15 mag, equal to the median magni- 
tude uncertainty in our real data. Simulations with scatter 
based on more realistic, variable measurement uncertaini- 
ties, including using the distribution of uncertainties from 
the real data, produce very similar results, with a slight 
broadening of the TFR and corresponding extension of the 
AMb TF versus log V TO t correlation, due to the tail of points 
with large uncertainties. 

We apply the scatter due to measurement errors after 
the magnitude cut, as the magnitude limit is approximately 
due to our inability to measure the properties of galaxies 
with real magnitudes fainter than the limit, rather than a 
direct cut based on our observed apparent magnitudes. 

'True' log Kot are generated for the simulated galaxies, 
using the 'true' Mb and the example TFR used through- 
out these simulations. These are converted into 'observed' 
log Kot by adding random offsets to mimic intrinsic scat- 
ter and measurement errors, both drawn from Gaussian dis- 
tributions. The measurement error standard deviation used 
was the median log V I0 t uncertainty in our real data, 0.04 
dex. The intrinsic scatter standard deviation used was 0.12 
dex, corresponding to the 0.9 mag of intrinsic scatter deter- 
mined in fitting our real field TFR (see Eqn.0. 

The above procedure has been repeated to generate 10 8 
simulated points. The resulting distribution of these simu- 
lated galaxies in the TFR plot is shown by panel (a) of Fig. [5] 
The distribution's properties are broadly similar to those of 
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the simpler examples in Fig.|H|(c) and (d), but with a much 
more realistic appearance. Note that the sharp cut-off at the 
bright end of the luminosity function produces a fairly well 
defined edge to the top of the TFR, while the apparent mag- 
nitude limit causes a gentler drop-off at the bottom of the 
TFR. This is because the wide redshift range of the sample 
smears out this limit in terms of absolute magnitude. 

Residuals of the simulated data points from the fiducial 
TFR are calculated in the same manner as for the real data. 
Panel (b) of Fig. |3] shows the distribution of the simulated 
galaxies in the AMb TF - log Kot plane. Note the clear cor- 
relation in this plot, despite the fact that the points were 
generated from a TFR parallel to the fiducial PT92 TFR, 
and with only symmetrical scatter in Mb and log Kot ■ The 
well defined sloping top-right edge of this AMb TF - log Kot 
distribution is a consequence of the bright cut-off of the un- 
derlying luminosity function. Slightly less well-defined, but 
still clear, is the sloped edge defining the bottom-left limit 
of the distribution. This is due to the fairly gentle cut-off 
at the faint end of the TFR, itself a result of the apparent 
magnitude limit. As this cut-off is gentler than at the bright 
end of the TFR, the slope in the AM S TF - log Kot plot is 
flatter, and thus the distribution broadens slightly towards 
low log V rot ■ 

The AMs TF - log Kot plot of our real data is shown in 
panel (c) of Fig. |5] These points display a correlation very 
similar to that seen in the simulated data sets. Note the sim- 
ilarity with the detailed simulation in panel (b). The only 
significant difference could be argued to be the existence of a 
few points at lower log Kot and higher AMs TF than the sim- 
ulated distribution. However, remember that the error dis- 
tribution in our simulation assumes a constant measurement 
uncertainty, whereas the real data shows a wide variation. 
Including a fraction of points with significantly larger mea- 
surement uncertainties causes the AMb TF - log Kot distri- 
bution to extend futher, to lower log Kot and higher AMs TF . 

The observed correlation of our data in the plot of 
AMb TF versus log Kot can therefore be easily explained as 
a result of the intrinsic correlation between the two axes, 
combined with scatter in the TFR and restrictions on the 
magnitude range, particularly due to the underlying lumi- 
nosity function. It is therefore clear that a correlation be- 
tween AMb TF and Kot does not necessarily imply a physi- 
cal effect, and in our case it appears to be entirely consistent 
with no Kot-dependent change in the TFR with redshift. 



with atotai = 1-2 mag and Omt = 1.0 mag, while for their 
full sample we obtain 

M B = (-6.4 ± 0.7) ■ log Kot + (-6.7 ± 2.3) . (12) 

with atotai = 1.0 mag and <7i nt = 0.8 mag. Our Tully-Fisher 
relations are thus in reasonable agreement. This is surpris- 
ing since we see no evidence for a flattening of the TFR 
slope, particularly when considering only the 'high-quality' 
data, while B04 claim to see a strikingly shallow slope. This 
may be in some part due to their interpretation of the cor- 
relation between AMb TF and Kot, which we have shown is 
potentially entirely due to correlated errors, rather than any 
underlying physical effect. They also use a different method 
of fitting the TFR, a bisector fit. This may be more easily 
biased to flatter slopes, by the effect of TFR scatter com- 
bined with restrictions on the magnitude range, than our 
purely inverse fit (see discussion in Section f3.3> . 

Apart from the disagreement over the evolution of the 
TFR slope, our results are very much in accord with those 
of B04. In particular we find very similar modest evolutions 
in the TFR offset with re dshift, of ~ 1 mag b y 2 = 1, also 
supported by the study of lBarden et alJ <l2003l) . 

Several other earlier studies of the Tully-Fisher re- 
lation at intermediate redshifts (e.g., iRix et all Il997l : 
ISimard fc Pritchetll99l) found significant luminosity evolu- 
tion with redshift, in agreement with ou r findi ngs. However, 
iDalcanton. Spergel fc Summersl (Il997l> and ISimard et alJ 
( 1999) argue that much, if not all, of the detected evolution 
is due to surface brightness selection effects. This is in agree- 
ment with the wo rk of Vogt and her co-workers IVoet et alJ 
1997, 2002) who find little or, more recently, no sig- 
nificant evolution in the TFR out to z — 1.3, once surface 
brightness selection effects have been accounted for. Further 
details of this latter study's selection and analysis methods 
will hopefully shed light on the difference between our re- 
sults and theirs. 

We will not attempt a more detailed analysis of the 
surface brightness selection effects in our study given the 
complexity of our selection. Nevertheless, for the above rea- 
sons and the arguments presented in Section ^. 21 we believe 
it is safe to take our measured luminosity evolution as an 
upper limit. In Section 2] we will use this limit to constrain 
the evolution of the star-formation rate of bright field spiral 
galaxies. 



3.4 Comparison with other studies 

The primary work with which we can compare is that of 
iBohm et alT(l2004l. B04). This is a similar study of field spi- 
rals between 0.1 < z < 1 with 36 'high quality' and 41 
'low quality' rotation velocity measurements. The internal 
extinction correction applied is the same as in this paper. 
In Fig. HOI we plot our data together with those of B04 4 . We 
also plot TFR fits to our data and their 'high quality' and 
full samples. Our fit to the B04 'high-quality' data is 

M B = (-9.5 ± 2.3) ■ log Kot + (0.2 ± 5.1) (11) 



4 obtained from the CDS catalogue archive 



4 SFR EVOLUTION WITH REDSHIFT 

We can consider what our best-fitting evolution in the TFR 
implies for the luminosity and star-formation evolution of 
field spirals. The implicit assumption in this section is that 
the £?-band luminosity of a spiral galaxies with a given rota- 
tion velocity evolves due to its star formation history only. 
For simplicity, we shall parameterise the luminosity and star- 
formation rate evolution as simple power-laws. 

Let us assume that the rest-frame B-band luminosity 
of field spirals evolves as a power law in (1 + z). Then, 

L B {z)=L B {0)-(l + zf (13) 

AMb(z) = Mb(z) - M s (0) = -2.5/3 log(l + z) (14) 

In Section |XU we find a best-fitting evolution with slope 
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Figure 10. A comparison of our TFR, with that of B04. The 
circular points and solid line are the same is in Fig. |S| as are 
the thin lines indicating the fiducial PT92 TFR. In addition the 
'high quality' and 'low quality' points from B04 are marked by 
large and small triangles respectively. The dashed line is a fit 
to the 'high quality' points, and the dotted line is a fit to the 
full B04 sample. Both of these fits are performed using the same 
algorithm used to fit our data, as described in Section 13, II The 
upper set of error bars in the top left corner indicate the 10th-, 
50th- and 90th-percentile uncertainties on the B04 'high quality' 
points, while the lower set give the same for our full TFR data 
set. 



AMb(z) w (-1.0 ± 0.5) • z mag 



(15) 



which implies Ait = 1.3 ± 0.6. Fitting directly with respect 
to log(l + z) gives the same value for /3. 

Before we proceed to compare with stellar evolution 
models, we need a way to relate (1 + z) to time. For the 
concordance cosmology, 



(1 + *) 



(16) 



is an excellent approximation for z < 1 and sufficiently ac- 
curate for our purposes up to z ~ 5. For a galaxy of age t age , 
which formed at tf we thus have 

{1 + z) ^( t l+^l)- - 813 (17) 



Using Eqn. 1171 to convert (1 + z) to t agc for vari- 
ous tf we can determine the expected value of (3 for any 
star-formation history usi ng the stellar evolution models of 
iBruzual fc Charlorl l|2003n . We model the SFR as 



SFRx 





(1 + *) 



' tf + tagc ' 
*0 / 



Z > Zf 
Z < Zf 



(18) 



For a galaxy with constant SFR (a = 0) and zt — 2 we find 
(3 a= o — —0.27. This corresponds to a brightening of 0.2 mag 
between z — 1 and today. A galaxy with constant SFR is 
therefore fainter in the past, opposite to what we measure. 
We therefore require a model with a SFR which declines 
with time, i.e., a > 



The difference in the power-law index between our best- 
fitting evolution and the constant SFR models is 



A/3 a=0 = At - A*=o = 1-6 ± 0.6 



(19) 



Choosing Zf = 5 only decreases this slightly to AA*=o = 
1.5 ± 0.6. We therefore work with Zf = 2 hereafter. 

In order to have some idea of of what value of a cor- 
responds to what is observed, we can assume that B-band 
luminosity is proportional to the current SFR, 

L B SFR 



Jfl,a=0 



SFR 



which implies a = AA*=o = 1.6 ± 0.6. 
form of for the SFR (i .e., SFR oc (1 + z, 
IBruzual fc Charlod i2003l) code gives 



(20) 



Putting this 
18 ) into the 



A*=i 



1.2 



A/3„=i. a = 0.1 



(21) 



nearly but not quite zero, as Eqn. l20l is not a perfect approxi- 
mation. However, we can improve our estimate by a straight- 
forward parameter search to determine the value of a which 
produces our observed Ait- We thus find a — 1.7 ± 1.1. 

In summary, we find that SFR(z) oc (1+z) 1 1 ' 1 in our 
sample of relatively bright (Mb > Mb) galaxies. As we ar- 
gued earlier, the luminosity evolution determined from our 
TFR evolution is likely to be an overestimate. The corre- 
sponding SFR evolution is therefore probably also an upper 
limit, and the evolution of the average SFR in luminous spi- 
rals may be slower than that found here. In contrast, stud- 
ies of the SFR den sity of the Universe at low redshift (e.g., 
iGalleeo etaITl 995) and as a function of look-back time (see, 
e.g.. lHeavens et all2 004 for a recent summary) generally in- 
dicate a very strong evolution up to z ~ 1. The precise rate 
of this evolution is somewhat controversial, but most studies 
agree that, pa rameterizing as SFR oc (1 + z) a , a ~ 3-4 (see 
iHopkinj 12001 for a compilation of SFR density evolution 
measurements) . 

This disparity suggests that the rapid evolution in the 
SFR density of the universe, observed since z ~ 1, is not 
driven by the evolution of the SFR in individual bright spi- 
ral galaxies, like those in our sample. However, it should be 
noted that, given the size of our sample and the intrinsic 
scatter in the TFR, the derived constraints on the evolu- 
tion of SFR in the bright spiral galaxy population are rela- 
tively weak. Nevertheless, this kind of approach seems very 
promising for future studies with samples of several hundred 
or even thousands of galaxies to z ~ 1. 



5 CONCLUSIONS 

We have measured confident rotation velocities for 89 field 
spirals with 0.1 < z < 1, and used these to examine the evo- 
lution of the TFR for these galaxies. The best-fitting TFR 
for the full sample has a slope which is entirely consistent 
with that found locally by PT92. There is an intercept offset 
of ~ 1 mag, such that our galaxies are brighter for a given 
Kot, but it is likely that this is a least partly due to the 
differing methods employed to measure the magnitudes and 
rotation velocities. We therefore only evaluate evolution in 
the TFR by making internal comparisons of our data. 

Fitting sub-samples binned by redshift indicates that 
the TFR intercept evolves by ~ 1 mag between z = 1 and 
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today, in the sense that more distant galaxies are brighter for 
a given V ro t- Plotting the residuals of our data from the local 
TFR against redshift confirms this trend in TFR intercept 
- higher-redshift galaxies are offset to brighter magnitudes. 
Fitting this correlation we find an evolution of — f .0 ± 0.5 
mag by z = 1, which we argue is an upper limit due to the 
selection effects present in our sample. 

We find no significant evidence for a change in TFR 
slope. Previous studies have used an observed correlation 
between the TFR residuals and Kot to argue that low mass 
galaxies have evolved significantly more than those with 
higher mass. However, we have demonstrated that such a 
correlation may be due solely to an intrinsic coupling be- 
tween Vrot scatter and TFR residuals, and thus does not 
necessarily indicate a physical difference in the evolution of 
galaxies with different Kot- 

Finally, we have used the stellar population models of 
iBruzual fc Charlod l)200ot) to interpret our observed TFR 
luminosity evolution in terms of a star-formation rate evo- 
lution. If the SFR in spiral galaxies had remained constant 
since z ~ 1, their B-band luminosity should have increased 
with time. However, we find the opposite trend, indicating 
that the SFR in these galaxies was larger in the past. We 
estimate SFR(z) oc (1 + z) 1 ' 7 1 for our sample. We argue 
that, given the likely selection effects, this is an upper limit 
on the SFR increase with look-back time. Our results there- 
fore suggest that the rapid evolution in the SFR density of 
the universe observed since z ~ I is not driven by evolution 
of the SFR in individual bright spiral galaxies. 

Even though we cannot yet place strong constraints on 
the evolution of the SFR in the bright spiral galaxy popu- 
lation, due to our relatively small sample and the intrinsic 
scatter in the TFR, our approach could be successfully ap- 
plied in ongoing and future surveys. A study with similar 
quality data for a sample of ~ 600 galaxies would be needed 
to reduce the error in a to ±0.4, providing a clear (~ 3<r) 
rejection of the hypothesis that the SFR density of the uni- 
verse and the SFR of the average spiral galaxy evolve at the 
same rate. 
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